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The Kibble-Zurek mechanism describes the evolution of topological defect structures like domain 
walls, strings, and monopoles when a system is driven through a second order phase transition. 
The model is used on very different scales like the Higgs field in the early universe or quantum 
fluids in condensed matter systems. A defect structure naturally arises during cooling if separated 
regions are too far apart to ‘communicate’ (e.g. about their orientation or phase) due to finite signal 
velocity. This results in separated domains with different (degenerated) locally broken symmetry. 
Within this picture we investigate the non-equilibrium dynamics in a condensed matter analogue, 
a two-dimensional ensemble of colloidal particles. In equilibrium it obeys the so called Kosterlitz- 
Thouless-Halperin-Nelson-Young (KTHNY) melting scenario with continuous (second-order like) 
phase transitions. The ensemble is exposed to a set of finite cooling rates covering roughly three 
orders of magnitude. Along this process, we analyze the defect and domain structure quantitatively 
via video microscopy and determine the scaling of the corresponding length scales as a function of 
the cooling rate. We indeed observe the scaling predicted by the Kibble-Zurek mechanism for the 
KTHNY universality class. 
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In the formalism of gauge theory with spontaneously 
broken symmetry, T.W.B. Kibble and colleagues postu¬ 
lated a cosmological phase transition during the cool¬ 
ing down of the early universe. This transition leads 
to degenerated states of vacua below a critical tempera¬ 
ture, separated or dispersed by defect structures as do¬ 
main walls, strings or monopoles m- In the course 
of the transition, the vacuum can be described via an 
TV-component, scalar order parameter (known as the 
Higgs field) underlying an effective potential 

V = a(j> 2 + b (<)> 2 — t]q) 2 (1) 

where a is temperature dependent, 6 is a constant and ijo 
is the modulus of (</>) at T = 0. For high temperatures, 
V has a single minimum at (f> = 0 (high symmetry) but 
develops a minimum ‘landscape’ of degenerated vacua be¬ 
low a critical temperature T c (e.g. the so called sombrero 
shape for N = 2). Cooling down from the high symme¬ 
try phase, the system undergoes a phase transition at T c 
into an ordered (low symmetry) phase with non-zero (6). 
For T <T C it holds 

<0) 2 = % 2 (l - T 2 /T 2 ) = v 2 (T) . (2) 

Caused by thermal fluctuations, one can expect that be¬ 
low T c , (</>) takes different non-zero values in regions 
which are not connected by causality. The question now 
arising concerns the determination of the typical length 
scale of these regions and their separation. For a finite 
cooling rate, fa is limited by the speed of propagating in¬ 
formation which is given by the finite speed of light defin¬ 
ing an ultimate event horizon. Independent of the nature 
of the limiting causality, Kibble argued that as long as the 
difference in free energy A F (of a certain system volume) 
between its high symmetry state (</>) = 0 and a possible 


finite value of {</>) just below T c is less than fcgT, the vol¬ 
ume can jump between both phases. The temperature 
at which A F = fcgT is called the Ginzburg temperature 
Tq and the length scale £d of the initial (proto)domains 
is supposed to be equal to the correlation length at that 
temperature: = £(Ig) [ 2]. 

The geometry of the defect network that separates the 
uncorrelated domains is given by the topology of the 
manifold of degenerated states that can exist in the low 
symmetry phase. Thus, it depends strongly on the di¬ 
mensionality of the system D and on the dimension N of 
the order parameter itself. Regarding the square root of 
Eqn. [2j the expectation value of a one-component order 
parameter (N = 1) can only take two different low sym¬ 
metry values ((/>) = ±i](T) (like e.g. the magnetization 
in a 2D or 3D Ising model): The manifold of the possible 
states is ‘disconnected’. This has a crucial effect if one 
considers a mesh of symmetry broken domains where (</>) 
is chosen randomly either +r] or —77. If two neighboring 
(but uncorrelated) domains have the same expectation 
value of {(j>), they can merge. Contrarily, domains with 
an opposite expectation value will be separated by a do¬ 
main wall in 3D (or a ‘domain line’ in 2D). At its center, 
the domain wall attains a value of (<j>) = 0, providing 
a continuous crossover of the expectation value between 
the domains (Fig. 00- Consider now N = 2: ((f>) can take 
any value on a circle, e.g. (c/i>) 2 = ( <j > x ) 2 + (4> y ) 2 = t? 2 (T) 
(all the order parameter values which are lying on the 
minimum-circle of the ‘sombrero’ are degenerated). Since 
the manifold of possible low symmetry states is now con¬ 
nected, {</>) can vary smoothly along a path (Fig. [TJd) . 
In a network of symmetry broken domains in two di¬ 
mensions at least three domains (in FigJTJo separated by 
dashed lines) meet at a mutual edge. On a closed path 
around the edge, the expectation value (</>} might be ei- 
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ther constant along the path (for a global, uniform 0) 
but can also vary by a multiple of 2ir (in analogy to the 
winding numbers in liquid crystals). In the first case, the 
closed path can be reduced to a point with (0) ^ 0 and 
no defect is built. If the path is shrunk in the second 
case, the field eventually has to attain (0) = 0 within the 
path and one remains with a monopole for D = 2 or a 
string for D = 3 [2j[3]. A condensed matter analogue 
would be a vortex of normal fluid with quantized circu¬ 
lation in superfluid helium. For N = 3 and D = 3, four 
domains can meet at a mutual point and the degenerated 
solutions of the low temperature phase lie on a sphere: 
(eft) 2 = (4>x) 2 + {(fry) 2 + (0z) 2 - If the held now again varies 
circularly on a ‘spherical path’ (all field arrows point ra¬ 
dially outwards), a shrinking of this sphere leads to a 
monopole in three dimensions Ell- 

W.H. Zurek extended Kibble’s predictions and trans¬ 
ferred his considerations to quantum condensed matter 
systems. He suggested that 4 He should intrinsically de¬ 
velop a defect structure when quenched from the normal 
to the superfluid phase mm- For superfluid 4 He, the 
order parameter 0 = |0| exp(?0) is complex with two 
independent components: magnitude |0| and phase 0 
(the superfluid density is given by |0| 2 ). A nontrivial, 
static solution of the equation of state with a Ginzburg- 
Landau potential yields 0 = 0o(r) exp {imp) where r and 
^ are cylindrical coordinates, n £ Z and "00(0) = 0. This 
solution is called a vortex line, topologically equivalent 
to a string for the case N = 2 we have discussed be¬ 
fore. In the vicinity of the critical temperature during a 
quench from the normal fluid to the superfluid state, 0 
will be chosen randomly in uncorrelated regions leading 
to a string network of normal fluid vortices. In condensed 
matter systems, the role of the limiting speed of light is 
taken by the sound velocity (in 4 He, the second sound). 
This leads to a finite speed of the propagation of order 
parameter fluctuations and sets a ‘sonic horizon’. 

Zurek argued that the correlation length is ‘frozen-out’ 
close to the transition point or even far before depend¬ 
ing on the cooling rate [4, 5 . Consider the divergence 
of the correlation length £ for a second-order transition, 
e.g. £ = £o | £ P' y where e = (T — T c )/T c is the reduced 
temperature. If the cooling is infinitely slow, the system 
behaves as in equilibrium: £ will diverge close to the tran¬ 
sition and the system is a mono-domain. For an instan¬ 
taneous quench, the system has minimal time to adapt 
to its surrounding: £ will be frozen-out at the beginning 
of the quench. For second order phase transitions, the 
divergence of correlation lengths is accompanied by the 
divergence of the correlation time r = tq |e| M which is 
due to the critical slowing down of order parameter fluc¬ 
tuations. If the time t it takes to reach T c for a given 
cooling rate is larger than the correlation time, the sys¬ 
tem stays in equilibrium and the dynamic is adiabatic. 
Nonetheless, for every finite but nonzero cooling rate, t 
eventually becomes smaller than r and the system falls 


out of equilibrium before T c is reached. This is the so 
called freeze-out time t. given when the correlation time 
equals the time it takes to reach T 0 : 

t = t(* ) . (3) 

The frozen out correlation length £ is then set at the 
temperature e of the corresponding freeze-out time: £ = 
£(e) = £(t). For a linear temperature quench 

e = (T- T c )/T c = t/r q (4) 

with the quench time scale r q , one observes t = 
( T o T q ) 1/(1+Al) and 

i=at) = Co(r q /r o r^ . (5) 

For the GL model (v = 1/2, /r = 1) one finds the scaling 
£ ~ T^ 4 while a renormalization group correction (v = 
2/3) leads to £ ~ t ^ 3 (4| [5]. 

A frequently used approximation is that when the adi¬ 
abatic regime ends at t before the transition, the corre¬ 
lation length can not follow the critical behavior until 
r again exceeds the time t when T c is passed. Given a 
symmetric divergence of r around T c , this is the time t 
after the transition. The period in between is known as 
the impulse regime in which the correlation length is as¬ 
sumed not to evolve further. A recent analytical investi¬ 
gation, however, suggests that in this period, the system 
falls into a regime of critical coarse graining [5J. There, 
the typical length scale of correlated domains continues 
to grow since local fluctuations are still allowed and the 
system is out of equilibrium. On the other side, numeri¬ 
cal studies in which dissipative contributions and cooling 
rates were alternatively varied before and after the tran¬ 
sition indicate that the final length scale of the defect and 
domain network is entirely determined after the transi¬ 
tion [7]. Several efforts have been made to provide exper¬ 
imental verification of the Kibble-Zurek mechanism in a 



FIG. 1. Emergence of defects in the Higgs-field which is il¬ 
lustrated with red vectors (shown in 2D for simplicity), a: 
For N = 1 and D = 3, domain walls can appear (strings for 
D = 2). b: For N = 2 and D = 3, nontrivial topologies 
are strings (monopoles for D = 2). The defects are regions 
where the order parameter 0 retains the high symmetry phase 
((0) = 0) to ’moderate’ between different degenerated orien¬ 
tations of the symmetry broken held. 
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variety of systems, e.g. in liquid crystals [3] (the transi¬ 
tion is weakly first order but the defect network can easily 
observed with cross polarization microscopy), superfluid 
3 He [H, superconducting systems uni , convective, intrin¬ 
sically out of equilibrium systems EL multiferroics El, 
quantum systems El, ion crystals El E, and Bose- 
Einstein condensates PH (the latter two systems contain 
the effect of inhomogeneities due to e.g. temperature 
gradients). A detailed review concerning the significance 
and limitations of these experiments can be found in E- 

In this experimental study, we test the valid¬ 
ity and applicability of the Kibble-Zurek mecha¬ 
nism in a two-dimensional colloidal model system 
whose equilibrium thermodynamics follow the microscop¬ 
ically motivated Kosterlitz-Thouless-Halperin-Nelson- 
Young (KTHNY) theory. This theory predicts a con¬ 
tinuous, two-step melting behavior whose dynamics, 
however, are quantitatively different from phenomeno¬ 
logical second-order phase transitions described by the 
Ginzburg-Landau model. We applied cooling rates over 
roughly three orders of magnitude for which we changed 
the control parameter with high resolution and homo¬ 
geneously throughout the sample without temperature 
gradients. Single particle resolution provides a quantita¬ 
tive determination of defect and domain structures dur¬ 
ing the entire quench procedure, and the precise knowl¬ 
edge of the equilibrium dynamics allows to determine 
the scaling behavior of corresponding length scales at the 
freeze-out times. In the following, we validate that the 
Kibble-Zurek mechanism can be successfully applied to 
the KTHNY universality class. 

DEFECTS AND SYMMETRY BREAKING IN 2D 
CRYSTALLIZATION 

The closed packed crystalline structure in two dimen¬ 
sions is a hexagonal crystal with sixfold symmetry. The 
thermodynamics of such a crystal can analytically be de¬ 
scribed via the KTHNY theory, a microscopic, two-step 
melting scenario (including two continuous transitions) 


FIG. 2. Sketch of a fivefold oriented (a) and sevenfold ori¬ 
ented (b) disclination. The red arrows illustrate the change 
in bond angle (blue) when circling on an anti clock-wise path 
around the defect. 


which is based on elasticity theory and a renormaliza¬ 
tion group analysis of topological defects [T8142CT] . In 
the KTHNY formalism, the orientationally long-range- 
ordered crystalline phase melts at a temperature T m via 
the dissociation of pairs of dislocations into an hexatic 
fluid, which is unknown in 3D systems. This fluid is char¬ 
acterized by quasi-long-range orientational but short- 
range translational order. In a triangular lattice, disloca¬ 
tions are point defects and consist of two neighboring par¬ 
ticles with five and seven nearest neighbors, respectively, 
surrounded by sixfold coordinated particles. At a higher 
temperature T), dislocations start to unbind further into 
isolated disclinations (a disclination is a particle with five 
or seven nearest neighbors surrounded by sixfold coordi¬ 
nated particles) and the system enters an isotropic fluid 
with short-range orientational and translational order. A 
suitable orientational order parameter is the local bond 
order field ip 6 (fj,t) = nj 1 e l69jk ^ t ' > = \ip 6 (fj,t)\e iej ^ 
which is a complex number with magnitude \i/jQ(rj,t)\ 
and phase Qj(t) defined at the discrete particle posi¬ 
tions fj. Oj{t) is the average bond orientation for a 
specific particle. The fc-sum runs over all rij nearest 
neighbors of particle j, and 9j k is the angle of the fcth 
bond with respect to a certain reference axis. If par¬ 
ticle j is perfectly sixfold coordinated (e.g. all 9j k (t) 
equal an ascending multiple of 7r/3), the local bond or¬ 
der parameter attains |'06(f = j,i)| = 1- A five- or sev¬ 
enfold coordinated particle yields > 0. The 

three different phases can be distinguished via the spa¬ 
tial correlation g 6 (r) = or temporal cor¬ 

relation g e (t) = ( , i/’g(0)V’6(t)) of the local bond order pa¬ 
rameter. For large r and t , ge(r,t) attains a finite value 
in the (mono)crystalline phase, decays algebraically in 
the hexatic fluid, and exponentially ~ exp(— r/£g) and 
~ exp (—t/re) in the isotropic fluid (US ?j . Unlike 
second-order phase transitions where correlations typi¬ 
cally diverge algebraically, the orientational correlation 
length ^6 and time tq diverge in the KTHNY formalism 
exponentially at T: 

~ exp(a |ep 1/2 ) and r 6 ~ exp(& |e| _1/2 ) , (6) 

where e = (T — T))/!), and a and b are constants [201122] . 
This peculiarity is the reason why KTHNY-melting is 
named continuous instead of second-order. In equilib¬ 
rium, the KTHNY scenario has been verified successfully 
for our colloidal system in various experimental stud¬ 
ies [23H25]. 

To transfer this structural 2D phase behavior into the 
framework of the Kibble-Zurek mechanism, we start in 
the high temperature phase (isotropic fluid) and describe 
the symmetry breaking with the spatial distribution of 
the bond order parameter. Since in 2D the local symme¬ 
try is sixfold in the crystal and the fluid, the isotropic 
phase is a mixture of sixfold and equally numbered five- 
and sevenfold particles (other coordination numbers are 
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FIG. 3. Orientational correlation time tq (experimental data 
and fit according to Eqn. [6j and the time t left until the 
transition temperature is reached for different cooling rates 
(colored straight lines, Eqn. [9j) as function of inverse temper¬ 
ature r (small r correspond to large temperatures and vice 
versa). The intersections define the freeze-out interaction val¬ 
ues f = r(t). 


extremely rare and can be neglected). During cooling, 
isolated disclinations combine to dislocations which for 
‘infinite’ slow cooling rates can annihilate into sixfold 
particles with an uniform director field. This is given by 
a global phase, characterizing the orientation of the crys¬ 
tal axis. Spontaneous symmetry breaking implies that all 
possible global crystal orientations are degenerated and 
the Kibble-Zurek mechanism predicts that in the pres¬ 
ence of critical fluctuations the system can not gain a 
global phase at finite cooling rates: Locally, symmetry 
broken domains will emerge which will have different ori¬ 
entations in causally separated regions. The final state 
is a polycrystalline network with frozen-in defects. As 
in the case of superfluid 4 He, ipe(fj,t) is complex with 
two independent components (N = 2). Consequently, 
we expect to observe monopoles in two dimensions. The 
phase of ipQ (r), t) is invariant under a change in the par¬ 
ticular bond angles of A 6jk{t) = ±mr/3 (n £ N) which 
is caused by the sixfold orientation of the triangular lat¬ 
tice. Similar to the Higgs field or the superfluid, one can 
not consider a closed (discrete) path in on which 

9j k (t) changes by an amount of ±7t/3 leaving the orien¬ 
tational field invariant. Reducing this path to a point, 
ipe (r, t) must tend towards zero at the center to maintain 
continuity. Since the orientational field is defined at dis¬ 
crete positions, the defect is a single particle marked as a 
monopole of the high symmetry phase. In fact, this coin¬ 
cides with the definition of disclinations in the KTHNY 
formalism [201 : The particle at the center is an isolated 
five- or sevenfold coordinated site. Fig. [2] illustrates this 
for a bond on a closed path. Going counter clock-wise 
around the defect, the bond angle changes by an amount 
of +7t/3 for a fivefold (Fig. io and by —n/3 for a sev¬ 
enfold site (Fig. §>)■ (In principle, also larger changes in 
9jk(t) are possible, e.g. for n = 2, a four- or eightfold ori¬ 


ented site, but these are extremely rare.) In KTHNY the¬ 
ory the monopoles (disclinations) combine to dipoles (dis¬ 
locations) which can only annihilate completely if their 
orientation is exactly antiparallel. At finite cooling rates 
they arrange in chains, separating symmetry broken do¬ 
mains of different orientation: Chains of dislocations can 
be regarded as strings or ‘2D domain walls’. 


COLLOIDAL MONOLAYER AND COOLING 
PROCEDURE 

Our colloidal model system consists of polystyrene 
beads with diameter a = 4.5 /im, dispersed in water and 
sterically stabilized with the soap sodium dodecyl sulfate 
(SDS). The beads are doped with iron oxide nanoparti¬ 
cles which results in a superparamagnetic behavior and 
a mass density of 1.7 kg/dm 3 . The colloidal suspension 
is sealed within a millimeter sized glass cell where sedi¬ 
mentation leads to the formation of a monolayer of beads 
on the bottom glass plate. The whole layer consists of 
> 10 5 particles and in a 1158 x 865 /jm 2 sub window 
ss 5700 particles are tracked with a spatial resolution 
of sub-micrometers and a time resolution in the order of 
seconds. The system is kept at room temperature and ex¬ 
empt from density gradients due to a months-long precise 
control of the horizontal inclination down to //rad. The 
potential energy can be tuned by an external magnetic 
field H applied perpendicular to the monolayer which in¬ 
duces a repulsive dipole-dipole interaction between the 
particles. The ratio between potential energy E mag and 
thermal energy k B T, 

r = Lmag = /ip (t m ) 3/2 (xH ) 2 

k B T 4nk B T ’ 1 ’ 

acts as inverse temperature (or dimensional pressure for 
fixed volume and particle number), n = 1/(Zq is the 2D 
particle density with a mean particle distance ao ~ 13 
H m, and x = 1-9 ■ 10 -11 Am 2 /T is the magnetic sus¬ 
ceptibility of the beads. F is the thermodynamic con¬ 
trol parameter: A small magnetic field corresponds to 
a large temperature and vice versa. Measured values of 
the equilibrium melting temperatures are F m « 70.3 for 
the crystal/hexatic transition and Tj ss 67.3 for the hex- 
atic/isotropic transition [25,. The cooling procedure is 
the following: We equilibrate the system deep in the 
isotropic liquid at To ~ 25 and apply linear cooling 
rates I = AT/At deep into the crystalline phase up to 
Tend ~ 100 thenceforward we let the system equilibrate. 
We perform different rates, ranging over almost three 
decades from f = 0.000042 1/s up to T = 0.0326 1/s. 
The slowest cooling rate corresponds to a quench time 
of ~ 19 days and the fastest one to ~ 40 minutes. We 
would like to emphasize that with the given control pa¬ 
rameter there is no heat transport from the surface as 
in 3D bulk material. The lack of gradients rules out a 














5 


temperature gradient assisted annealing of defects which 
might be present in inhomogeneous systems. 


STRUCTURE AND DYNAMICS OF DEFECTS 
AND DOMAINS 


The key element of the Kibble-Zurek mechanism is a 
frozen-out correlation length £ as the system falls out of 
equilibrium at the freeze-out time t. For slow cooling 
rates, the system can follow adiabatically closer to the 
transition (large £) than for fast rates where the systems 
reaches the freeze-out time earlier (small £). To find t , we 
determine the orientational correlation time tq accord¬ 
ing the KTHNY theory by fitting g 6 (t) ~ exp(—f/r 6 ) in 
the isotropic fluid for independent equilibrium measure¬ 
ments. The data for as well as a fit with 

r 6 = r 0 exp(6 r |1/T - 1/T c p 1/2 ) (8) 


is shown on the right axis of Fig. [3] The time left to the 
isotropic-hexatic transition is given by 

t = (Fi - r)/f (9) 

and is also plotted in Fig. [3] (left axis) for various cool¬ 
ing rates including the slowest and the fastest one. The 
points of intersection 


t = tq exp 



i/r(i,r) - i/r c 



( 10 ) 



FIG. 4. Defect number density p and average domain size 
(A) (in units of Oq) as a function of T (oc 1/T) during cooling 
from small T (= ‘hot’ on the left side) to large F (= ‘cold’ on 
the right side). The curves cover the complete range of cooling 
rates from f = 0.000042 to f = 0.033 and are averaged within 
an interval Ar = 0.4. Big dots mark the freeze-out tempera¬ 
tures f = F(t) (colored correspondingly to f). Open symbols 
show the equilibrium melting behavior (lines are guide to the 
eye). 


define the freeze-out ‘temperatures’ f = F(f). 

The length scale of the defect network can be measured 
by the concentration p of defects (counting all not sixfold 
coordinated particles) in the ipe(f,t)-&e\ d. Fig. 0 (up¬ 
per plot) shows the evolution of p for the same cooling 
rates T as in Fig. pi as well as for the equilibrium (melt¬ 
ing) behavior [25^ One recognizes that the course of 
p deviates from the equilibrium case in advance of the 
isotropic/hexatic transition at ss 67.3. This happens 
at different times for distinct cooling rates and marks the 
end of the adiabatic regime. Within the noise, deviations 
from the equilibrium behavior start at the temperature T 
given by the freeze-out time t (big colored dots). Beyond 
the adiabatic regime the defect density decreases which 
is an indication of critical coarse graining as predicted 
in [SJ. At Ti (and also F m ) the slope of the curves in¬ 
creases with decreasing cooling rate indicating a further 
evolution but the system cannot perform critical fluctu¬ 
ations. 

The domain structure, on the other hand, can be 
characterized quantitatively by analyzing symmetry bro¬ 
ken domains with similar phase of ipe(rj,t). According 
to [26] . we define a particle to be part of a symmetry 
broken domain if the following three conditions are ful¬ 
filled for the particle itself and at least one nearest neigh¬ 
bor: 1) The magnitude \i/)6(rj : t)\ of the local bond order 
parameter must exceed 0.6 for both neighboring parti¬ 
cles, 2) the bond length deviation of neighboring par¬ 
ticles is less than 10% of the average particle distance 
ao, and 3) the variation in the average bond orienta¬ 
tion A &ij(t) = \Im[ip6(fi)] — Im[tpe(fj)]\ of neighbor¬ 
ing particles i and j must be less than 14° (less than 
14°/6 in real space). Simply connected domains of par¬ 
ticles which fulfill all three criteria are merged to a local 
symmetry broken domain. If a particle does not satisfy 
these conditions in respect to a neighboring particle, it 
is assigned to the high symmetry phase (almost all de¬ 
fects are identified as such due to their small value of 
Fi s-0 (lower plot) shows the evolution of 
the ensemble average domain size (A) as a function of 
r. We observe a behavior analogue to p\ Domain for¬ 
mation significantly deviates from the equilibrium case 
before F^, namely around the freeze-out temperature T 
of the corresponding cooling rates I\ To compare both 
networks in the following, we define the dimensionless 
lengths £def = p~ 1/2 and £ dom = ((-4) /ao) 1/2 which dis¬ 
play the characteristic length scales in units of ao- 

Colloidal ensembles offer the unique possibility to mon¬ 
itor the system and its domain and defect structure on 
single particle level. Fig. [5] illustrates both (left row for 
defects, right row for domains) at the freeze-out temper¬ 
ature T for the fastest (a,b) and the slowest (c,d) cooling 
rate. For T = 0.0326 1/s (a,b) where t is already reached 
at T = 30.3, the defect density is large as is the number of 
high symmetry particles. However, there is a significant 
number of sixfold coordinated particles and a few orienta- 
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FIG. 5. Snapshot sections of the colloidal ensemble 
(992 x 960 /im 2 ,« 4000 particles) illustrating the defect (a,c) 
and domain configurations (b,d) at the freeze out temperature 
f for the fastest (a,b: f = 0.0326 1/s, f ss 30.3) and slowest 
cooling rate (c,d: f = 0.000042 1/s, f rts 66.8). The defects 
are marked as follows: Particles with five nearest neighbors 
are colored red, seven nearest neighbors green and other de¬ 
fects blue. Sixfold coordinated particles are colored grey. Dif¬ 
ferent symmetry broken domains are colored individually and 
high symmetry particles are displayed by smaller circles. 


FIG. 6. Snapshot sections of the colloidal ensemble il¬ 
lustrating the defect (a,c) and domain configurations (b,d) 
after quasi-equilibration of the system for the fastest (a,b: 
F = 0.0326 1/s, r en d « 105) and slowest cooling rate (c,d: 
f = 0.000042 1/s, F cn d 98). The system size and the label¬ 
ing of defects and domains is the same as in Fig. [5] 

to be a monodomain but its observation lies beyond ex¬ 
perimental accessible times for our system. 


tionally ordered domains (to accord for finite size effects, 
we will exclude domains which hit the border of the field 
of view when evaluating £dom at T). At this point the 
length scales are £def = 1.56±0.01 and £dom = 1.56±0.03. 
For the slowest cooling rate F = 0.000042 1/s (c,d) where 
r = 66.7, the mean distance between defects as well as 
the typical domain size is significantly larger compared 
to the fastest cooling rate. We observe £def = 2.36±0.07, 
and £dom = 2.30 ± 0.09. 

To allow relaxation of the defect and domain structure 
after the freeze-out time we keep the temperature 
constant after T en d « 100 is reached. Fig. [6] shows the 
defect and domain configurations after an equilibration 
time of « 5 hours for the fastest cooling rate (a,b) where 
the quench time was S3 40 minutes, and after an equilibra¬ 
tion time of S3 3 days for the slowest rate (c,d) where the 
quench time was ss 19 days. The different length scale 
of the defect and symmetry broken domain network in 
respect to the cooling rate is clearly visible: While we 
observe a large number of domains for fast cooling, slow 
cooling results in merely two large domains separated by 
a single grain boundary. The final evolution will be given 
by classical coarse graining. The ground state is known 


SCALING BEHAVIOR 


The main prediction of the Kibble-Zurek mechanism 
is a power law dependence of the frozen-out correlation 
length £ as function of r q (Eqn. [ 5 ]) which results from the 
algebraic divergence of the correlation, presuming Eqn. [3] 
and Eqn.[4j In KTHNY-melting, £ 6 and r 6 diverge expo¬ 
nentially and one has to solve Eqn.[l0]to find the implicit 
dependency t(T). We did this numerically for discrete 
values in the range 4 • 10 -5 < F < 4 • 10 -2 and deter¬ 
mined the frozen-out orientational correlation length ^6 
for a scaling t s /tb = c with the dynamical exponent 
z J22]- Here, tb ~ 171.6 s is the Brownian time which is 
the time a single particle needs to diffuse its own diame¬ 
ter. Using Eqn. [i one finds with F(f, T) from Eqn. [ 9 ] the 
expression 


&(T) = 




r e - r, + ft(f) 
r c r, - r c ft(f) 


-1/2N 


( 11 ) 

This function is plotted in Fig.[7]for z = 4.5 and c = 0.83 
(red curves) on a double logarithmic scale together with 
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£def and £dom at the freeze-out temperature f. We 
find very good agreement. Nonetheless, we fit for com¬ 
parison the data via an algebraic scaling (blue dotted 
lines) of the form f(T) = aT~ K , for which we observe 
/vdef = 0.061 ± 0.001 and Kdom = 0.061 ± 0.002. The 
data are compatible with the algebraic decay only for in¬ 
termediate cooling rates. The deviations from standard 
Kibble-Zurek mechanism for systems with second order 
transitions are in line with the temperature-quenched 2D 
XY model m having also non algebraic divergences of 
the correlation length in equilibrium and being thus in a 
similar universality class. The small algebraic exponent 
k can be explained by the relatively large value of the dy¬ 
namical exponent z which regulates the ‘slope’ of ^(F) 
(in [22], a value z = 2.5 was proposed for the hard-disk 
system). This is due to quite long correlation times in 
this colloidal system (see Fig. [3j) which are caused by its 
overdamped dynamics. Note that the ‘sonic horizon’ is 
set by the sound velocity of the colloidal monolayer (and 
not the solvent) being ~mm/s, six orders of magnitude 
slower compared to atomic systems. 

CONCLUSIONS 

We presented a colloidal model system, where struc¬ 
ture formation in spontaneously symmetry broken sys¬ 
tems can be investigated with single particle resolution. 
The theoretical framework is given by the Kibble-Zurek 
mechanism which describes domain formation on differ- 



10 4 10 3 10 2 

f 


FIG. 7. The length scale of the defect £def and domain 
network £dom is plotted as a function of the cooling rate F 
(open symbols). Red lines are numerical solutions of the tran¬ 
scendental equation following the freeze-out condition for the 
KTHNY-like divergences (see text for definition). For com¬ 
parison, dashed blue lines are power law fits predicted by 
standard Kibble-Zurek mechanism which show the same al¬ 
gebraic exponent n « 0.06 for £def and (dem¬ 


ent scales like the Higgs field shortly after the Big Bang 
or the vortex network in 4 He quenched into the superfluid 
state. Along various cooling rates, we analyzed the devel¬ 
opment of defects and symmetry broken domains when 
the systems falls out of equilibrium and fluctuations of 
the order parameter can not follow adiabatically due to 
critical slowing down. While 2D melting in the colloidal 
monolayer is described by KTHNY theory where the di¬ 
vergence of the relevant correlation lengths in equilibrium 
is exponential (rather than algebraic as typically found in 
3D systems), the central idea of the Kibble-Zurek mech¬ 
anism still holds and the scaling of the observed domain 
network is correctly described. Implicitly this shows that 
existence of grain boundaries can not solely be used as 
criterion for first order phase transitions and nucleation 
or to falsify second/continuous order transitions since 
they naturally arise for non-zero cooling rates. Those 
will always be present on finite time scales in experiment 
and computer simulations after preparation of the sys¬ 
tem. 
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